Bose atoms in a trap: a variational Monte Carlo formulation for the universal 
behavior at the Van der Waals length scale 
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We present a variational Monte Carlo (VMC) formulation for the universal equations of state at 
the Van der Waals length scale [B. Gao, J. Phys. B 37, L227 (2004)] for N Bose atoms in a trap. The 
theory illustrates both how such equations of state can be computed exactly, and the existence and 
the importance of long-range atom-atom correlation under strong confinement. Explicit numerical 
results are presented for A'' = 3 and 5, and used to provide a quantitative understanding of the 
shape-dependent confinement correction that is important for few atoms under strong confinement. 

PACS numbers: 05.30. Jp,21.45.-f v,03.75.Nt 



INTRODUCTION 



Few atoms in a trap, which can in principle be real- 
ized, e.g., through a Mott transition from a degenerate 
quantum gas in a optical lattice P, Q j is a fundamental 
system for studying atomic interactions and correlations. 
It has many features of a many-atom system, yet still suf- 
ficiently simple to be amenable to a number of different 
theoretical approaches that can check against and learn 
from each other. These methods include, for example, 
the Monte Carlo methods (see, e.g^ II IIM IE P ), and 
the hyperspherical methods Iili[lS|n[ni|lllli. 

Unlike the problem of two identical atoms in a sym- 
metric harmonic trap, which can be reduced to a 1-D 
problem and be solved exactly [3 13 IHIlE IH , the 
problem of N atoms in a trap {N > 2) is much more 
complex. Partly due to this complexity, and partly due 
to our previously limited understanding of atomic inter- 
action based on the effect-range theory [3 1^ , existing 
studies of N atoms in a trap, with virtually no exception, 
have relied on potential models that do not reflect the 
true nature of atomic interaction except for the scatter- 
ing length. In other words, they provide understanding 
only at the level of shape-independent approximation. 
This approximation clearly breaks down for dense sys- 
tems with ~ 1 or greater |3|. Here p = N/V 
is the atomic number density, and Pe = (mC^I'fi'Y^^ is 
the length scale associated with the Van der Waals inter- 
action, —Cq/t^, between atoms. For an inhomogeneous 
system of atoms in a trap, the shape-independent approx- 
imation may also break down, even for fairly small values 
of p/3|, due to strong confinement 17, 18, 19, 24], as 
characterized by ag/aho ^ 1, where oq is the s wave scat- 
tering length, and aho = {h/muY^'^ is the length scale 
associated with the trapping potential. This effect, which 
we call the shape-dependent confinement correction [23] . 
can be understood qualitatively as due to the energy de- 
pendence of the scattering amplitude, which is always 



shape-dependent [11 13 111 E IH ■ 

Going beyond the shape-independent approximation 
requires understandings of atomic interaction and cor- 
relations at shorter length scales. Fortunately, uni- 
versal properties persist because atoms have the same 
types of long-range interactions, such as —Cq/t^ for 
atoms in ground state. The development of the angular- 
momentum-insensitive quantum defect theory (AQDT) 
pSL l26l I27L l28l ] has led both to a systematic understand- 
ing of atomic interaction of the types of — C„/r" with 
n > 3, and to a methodology for uncovering and study- 
ing universal properties at different length scales for two- 
atom, few- atom, and many-atom systems 22^.23.. 27. 29J. 
This work illustrates how this method can be imple- 
mented in a variational Monte Carlo (VMC) formulation 
that gives basically exact results for the iV-atom univer- 
sal equations of state at length scale (3q [H |2^ . Explicit 
numerical results are presented for three and five Bose 
atoms in a symmetric harmonic trap. They provide both 
samples of benchmark (basically exact) results for few 
atoms in a trap and a quantitative understanding of the 
shape-dependent confinement correction |24| . In the pro- 
cess of achieving these results, we also show that atoms 
in a trap have long-range correlation that becomes im- 
portant under strong confinement. 

Our VMC formulation for N Bose atoms in an external 
potential, which differs from existing formulations jj] in 
its choice of correlation function, is presented in Sec. ^ 
The universal equations of state at length scale (3q are 
discussed in Sec. IIIII with explicit numerical results for 
three and five Bose atoms in a symmetric harmonic trap 
presented in Sec. IIVI Conclusions are given in Sec. [V] 



II. VMC TREATMENT OF TV BOSE ATOMS IN 
AN EXTERNAL POTENTIAL 



The relative merit of different Monte Carlo methods 
are well documented We choose here the variational 
Monte Carlo method (VMC) for a number of reasons 

(a) VMC always works, for bosons, fermions, or excited 

states, provided one picks the right trial wave function. 
* Email |FgaoQphyslcs.utoIedo.edu| Homepage |http : //bgaowww . physics (Tlcb )3iV(M.(edp rovides the most transparent understanding 
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of many-body wave function, and is thus the best for 
conceptual purposes, (c) The advantages of other Monte 
Carlo methods Q , such that being a "black box" , mostly 
disappear when applied to fermions or to the excited 
states of a many-body system, (d) The result of VMC 
can always be used as the starting point upon which fur- 
ther adjustment or relaxation of wave function can be 
allowed, if at all necessary. More specifically, it can be 
used to fix the nodal structure and provide importance 
sampling 

The difficulty, or the challenge of VMC, is in choosing 
a proper trial wave function. Otherwise no converged 
result would be obtained, as reflected in the fact that 
the variance of energy would be of the same order of, 
or greater than, the average value being evaluated. The 
same challenge can, however, also be regarded as an op- 
portunity, as it forces one to understand the wave func- 
tion. 

Consider an A^-atom Bose system described by the 
Hamiltonian 
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(1) 



where Vext describes the external "trapping" potential, 
and v{r) represents the interaction between atoms that 
has a behavior of v{r) — > —C'e/r^ in the limit of large r. 

Such an iV-atom Bose system has of course many differ- 
ent states. We focus ourselves here on the lowest gaseous 
Bose-Einstein condensate (BEC) state, which can be de- 
fined as the state that evolves from the lowest TV-free- 
particle state in a trap as one turns on an atomic inter- 
action with positive scattering length. For this particular 
state, we take the variational trial wave function to be of 
Jastrow form 'sS] 



■ N 
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.i=l 



N 

n 



(2) 



It is straightforward to show that the expectation value 
of energy for such a state can be written as 

Z7 — J 



/ dr^*^£;L(ri,r2,r3) 



(3) 



where the integrations are over all iV-atom coordinates, 
and El is a local energy that can written as the sum of 
three terms whose contributions to the energy depend on 
the 1-body, 2-body, and three-body correlation functions, 
respectively: 



i?L = 4'^(ri) + i?r(ri,r2) 
Here 
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F(ri2)F(ri3) 



[ViF(ri2)] • [ViF(ri3)] . (9) 



Once (f) a-nd F are chosen, Eq. Q can be evaluated us- 
ing Metropolis Monte Carlo method (see, e.g., and 
the variational parameters are then varied to find the 
stationary energies. 

The success, or the failure, of a VMC calculation de- 
pends exclusively on the proper choice of the wave func- 
tion. The choice of is fairly standard and is based on 
the independent-particle solution in the external poten- 
tial. The choice of F is less obvious, and depends on 
the understanding of atom-atom correlation in a trap. 
Our choice of F is based the following physical consid- 
erations, (a) Atom-atom correlation at short distances 
is determined by two-body interaction, (b) Atoms in a 
trap can have long-range correlation that becomes im- 
portant under strong confinement, as suggested by our 
recent work on two atoms in a trap . Specifically, we 
choose our F as 



Aux{r)/r 



r < d 
r> d 



Here u{r) satisfies the Schrodinger equation: 



TO dr2 



+ v(r) - A 



u\{r) = , 



(10) 



(11) 



for r < d. 7 is the parameter characterizing the long- 
range correlation between atoms in a trap, with 7 = 
(meaning F = 1 for r > d) corresponding to no long- 
range correlation. Both d and 7 are taken to be varia- 
tional parameters, in addition to the variational param- 
eters associated with the description of 0. The parame- 
ters A and A are not independent and are determined by 
matching F and its derivative at d. 

The key difference between our choice of F and the 
standard choices 0, in addition to the systematic treat- 
ment of atomic interaction to be discussed in the next 
section, is the allowance for the long-range correlation 
characterized by parameter 7 [T9| . One can easily verify 
that regardless the model potential used for v (such as 
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the hard sphere potential), a choice of F without long- 
range correlation, such as 



F{r) = 1 - aoA 



(12) 



would not have led to converged VMC results under 
strong confinement. This explains why the existing 
Monte Carlo results for few atoms under strong confine- 
ment have come from diffusion Monte Carlo (DMC) Igl, 
but not from VMC, which was successful for weak con- 
finements 



III. UNIVERSAL EQUATION OF STATE AT 
THE VAN DER WAALS LENGTH SCALE FOR N 
BOSE ATOMS IN A SYMMETRIC HARMONIC 
TRAP 

For any state in which the atomic interaction at 
the average atomic separation is well represented by 
—Ce/r^, which for N Bose atoms in a trap implies 
p/3g ~ N(0(i/ah„)^ <^ 10, its energy follows a univer- 
sal behavior |23, 123 that is uniquely determined by the 
trapping and the Van der Waals potentials, independent 
of the interactions at short distances except through a 
parameter that can be taken either as the short range 
K matrix K'^ 25] or the s wave scattering length uq. 
Within the VMC formulation, this can be understood by 
noting the for such diffuse states, the solution ux{r) of 
Eq. ((TT)l . wherever it has an appreciable value j2^, is 
given by [H El El HI 



oirs) 



(13) 



Here B is a normalization constant 
universal AQDT reference functions for 



fc(6) 



are 



and g'^^i^ 
~Ce/r^ type of 



They depend on r only through 
r//36, and on energy only through 



potential |22l |2£ 
a scaled radius rg 

a scaled energy = X/se, where se = {U'^ /m){l/ (3q)'^ 
is the energy scale associated with the Van der Waals 
interaction. K'^ is a short-range K matrix 1251 that is 
related to the s wave scattering length oq bv |27ll2^ 



ao//36 



,26r(i-&) 



r(i + fe) 



X= + tan(7r6/2) 
X= - tan(7r6/2) ' 



(14) 



where h = l/(n — 2), with n ~ Q. Note that while K'^ 
and ao are related to each other, by propagating the wave 
function in the Van der Waals potential from small to 
large distances |23, E3| , they have considerably different 
physical meanings. K'^ is a short-range parameter that 
is directly related to the logarithmic derivative of the 
wave function coming out of the inner region, a re gion 
where atomic interaction may differ from —C^jr^ |25j |. 
flo is determined by the asymptotic behavior of the wave 
function at large distances. The universal behavior is 
conceptually easier to understand in terms of if^, as it 
simply implies that for any state in which the probability 
for finding particles in the inner region is small, the only 



role of the inner region is in determining the logarithmic 
derivative of the wave function coming out of it. Our 
results are presented in terms of a scaled ao parameter 
only to facilitate connections with existing models and 
understandings. 

When ua, as given by Eq. (|13|l . and therefore F, de- 
pend on the interactions of shorter range than /3g only 
through K'^ or a scaled ag, so do the overall wave func- 
tion and the energy of the A^-atom Bose system. For 
an inhomogeneous system of atoms in a trap, the energy 
depends of course also on the trap configuration. To be 
specific, we consider here atoms in a symmetric harmonic 
trap, characterized by 



(15) 



where w is the trap frequency. The corresponding 
independent-particle solution suggests 



r) = exp [-a{r/aho)'^] , 



(16) 



where a is chosen as one of the variation parameters, 
in addition to parameters d and 7 used to characterize 
the correlation function F. From this combination of (j) 
and F, the resulting VMC energy per particle, properly 
scaled, can be written as 



E/N 



l( 



(17) 



Where O is a universal function that is uniquely deter- 
mined by the number of particles, the exponent of the 
Van der Waals interaction {n — 6), and the exponent of 
the trapping potential (2 for the harmonic trap). The 
strengths of interactions, as characterized by Cg and uj, 
play a role only through scaling parameters such as f3e 
and aho- 

Equation 117|l . which is one example of what we call 
the universal equation of state at length scale /Jg, can 
also be defined, independent of the VMC formulation, 
using the method of effective potential as in Ref. . It 
is a method of renormalization in the coordinate space 
to eliminate all length scales shorter than (3^. The same 
procedure in VMC corresponds simply to using Eq. H13() 
for all r < d. The function fi, following this proce- 
dure, is rigorously defined for all values of ao/aho and 
for all f3e/aho > 0. An A^-atom Bose system in a sym- 
metric harmonic trap and in the lowest gaseous BEC 
state can be expected to follow this universal behavior 
for Pe/o-ho <'~ 2/iV^/^, beyond which the interactions of 
shorter range, such as —Cs/r^, can be expected to come 
into play. 

It is worth noting that the parameter Pe/aho in 
Eq. I|17|) plays a similar role, for atoms in a trap, as p(3q 



for homogeneous systems [22, 123 • The latter parame- 
ter is not used here obviously because p is not uniform, 
but its order of magnitude is still related to Pe/aho by 
p/3g ~ N{f3Q/ahoY ■ When either parameter goes to zero, 
the universal equations of state at length scale /3g can be 
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FIG. 1: The universal equation of state for three atoms in a 
symmetric harmonic trap as a function of ao/aho, compared 
to the DMC results of Blume and Greene for hard spheres 0- 



expected to go to the shape-independent results as ob- 
tained by Blume and Greene for particles in a tra p an d 
by Giorgini et al. for homogeneous systems p3. l23l |. 



IV. RESULTS FOR FEW BOSE ATOMS IN A 
SYMMETRIC HARMONIC TRAP 

The formulation in previous sections is applicable to 
any number of atoms. We present here explicit numerical 
results for few Bose atoms in a symmetric harmonic trap. 
This is not only because such calculations are less numer- 
ically intensive than for larger number of atoms, but also 
because before N gets sufficiently large that p(3^ ^ 1, the 
shape-dependent confinement correction is actually more 
important for smaller number of particles |24j . 

Figure ^ illustrates the equation of state for three 
atoms in a symmetric harmonic trap. It is a function 
of two variables that we plot here as a set of functions 
of ao/a^o for different values of /Je/oho- The results 
for Pe/aho — 0.001 show that, as expected, the uni- 
versal equation of state at length scale /Jg does eventu- 
ally approach a shape-independent result in the limit of 
Pe/oLho ^ 0, and are in excellent agreement with the 
DMC results of Blume and Greene for hard spheres 
The results for /^e/a^o = 0.01 and Pe/aho = 0.1 il- 
lustrate the shape-dependence of the equation of state 
due to the Van der Waals interaction. They show that 
even for relative small p/3g, which is of order of 10~^ for 
Pe/oLho — 0.01, the shape-dependent correction can be- 
come quite appreciable under strong confinement. This 
correction, which we call the shape-dependent confine- 
ment correction |24[, can be understood qualitatively as 
due to ener gy depe ndence of the two-body scattering am- 
plitude [T^ll7lll8LIT9ll23 |. which becomes significant for 
large scattering lengths. To put our results in perspec- 
tive, we note that a recent experiment on two atoms in a 




FIG. 2: The parameter 7, characterizing the long-range atom- 
atom correlation, for three atoms in a symmetric harmonic 
trap, as a function of ao/aho- 
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FIG. 3: The universal equation of state for five atoms in a 
symmetric harmonic trap as a function of a^/aho, compared 
to the DMC results of Blume and Greene for hard spheres _6j . 



symmetric harmonic trap is already exploring the region 
close to /Sg/a/io ~ 0.1 2]. 

Figure El shows that the parameter 7, characterizing 
the long-range atomic correlation, for three atoms in a 
symmetric harmonic trap. It is clear that 7 can become 
quite large under strong confinement, ao/aho 1- Not 
surprisingly, a variational wave function that does not 
incorporate this long-range correlation explicitly would 
fail under such conditions. 

Figure [S] shows the equation of states for five atoms 
in a symmetric harmonic trap. Compared to the re- 
sults for three atoms, the shape-dependent corrections 
can be seen to be less significant, confirming the conclu- 
sion from Ref. that the shape-dependent confinement 
correction is more important for smaller number of par- 
ticles than for larger number of particles. The long-range 
atom- atom correlation is again very important, as shown 
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nance [3j, |35| , the same results would apply, provided it 
is a broad Feshbach resonance with a width much greater 
than the energy scale, sr, associated with the Van der 
Waals potential [MIstIIM Isgf . 



CONCLUSIONS 



FIG. 4: The parameter 7, characterizing the long-range atom- 
atom correlation, for five atoms in a symmetric harmonic trap, 
as a function of ao/aho- 

TABLE I: Selected data of of energy per particle, in units of 
Tiu!, as a function of ao/aho for three atoms in a symmetric 
harmonic trap. The number in the parenthesis represents the 
variance in the last digit. 



ao/aho 


HS" 


Pe,/aho = 0.001 


Pe/aho = 0.01 


Pe/aho = 0.1 


0.433 


1.851 


1.851(2) 


1.911(2) 


1.957(1) 


0.866 


2.233 


2.237(1) 


2.327(1) 


2.411(1) 


1.732 


3.107 


3.110(1) 


3.235(1) 


3.375(1) 


2.598 


4.154 


4.162(1) 


4.301(2) 


4.426(1) 



"DMC results for hard spheres from Ref. 



in Fig. 21 

Some specific data points shown in Fig. [21 and 0] are 
tabulated in Tables HI and Ull for the convenience of future 
comparisons. They represent samples of basically exact 
results for few Bose atoms in a symmetric harmonic trap. 
The effects of atomic interactions with shorter ranges 
than Pq would come into play only for states with en- 
ergies that are much further away from the threshold (ei- 
ther below or above). If the scattering length parameter 
used here is achieved by tuning around a Feshbach reso- 



In conclusion, we have presented a VMC formulation 
for the universal equations of state at the length scale (3q 
for N Bose atoms in a symmetric harmonic trap. We have 
also shown that atoms under strong confinement have sig- 
nificant long-range correlation of the form of . Since an 
independent-particle model, such as the Hartree-Fock ap- 
proximation, corresponds to a variational method based 
on a wave function with F ^ 1, the fact that F, for 
atoms under strong confinement, deviates significantly 
from 1 everywhere implies that any independent-particle 
model is likely to fail for such systems. The results for 
TV = 3 and 5 provide a quantitative understanding of the 
shape-dependent confinement correction, which is impor- 
tant for a small number of particles under strong confine- 
ment [2i |. 

We are extending our calculations to larger number 
of particles to study universal behaviors, for both ho- 
mogeneous and inhomogeneous systems, in the region of 
p/?g ~ I, where shape-dependence is obviously impor- 
tant [22, 123 • We are also extending our methodology to 
other states of few-atom and many-atoms systems. They 
include not only the excited states with higher energies 
than the lowest gaseous BEC states, but also the liquid 
states with lower energies [23, . 



TABLE IL Selected data of of energy per particle, in units 
of hu, as a function of ao/aho for five atoms in a symmetric 
harmonic trap. The number in the parenthesis represents the 
variance in the last digit. 

ao/aho HS" Pe/aho = 0.001 ft/flho = 0.01 Pe/aho = 0.1 



0.433 2.115 


2.116(1) 


2.159(1) 


2.210(1) 


0.866 2.720 


2.722(1) 


2.791(1) 


2.844(1) 


1.732 4.018 


4.019(2) 


4.101(2) 


4.163(1) 


2.598 5.560 


5.561(1) 


5.657(2) 


5.729(1) 


"DMC results for hard spheres 


from Ref. |||. 
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